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ABSTRACT: The possibilty of performing high-rate calorimetry with a slow scintillator crystal is 
studied. In this experimental situation, to avoid pulse pile-up, it can be necessary to base the energy 
measurement on only a fraction of the emitted light, thus spoiling the energy resolution. This effect 
was experimentally studied with a BGO crystal and a photomultiplier followed by an integrator, 
by measuring the maximum amplitude of the signals. The experimental data show that the energy 
resolution is exclusively due to the statistical fluctuations of the number of photoelectrons con- 
tributing to the maximum amplitude. When such number is small its fluctuations are even smaller 
than those predicted by Poisson statistics. These results were confirmed by a Monte Carlo simu- 
lation which allows to estimate, in a general case, the energy resolution, given the total number of 
photoelectrons, the scintillation time and the integration time. 
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1. Introduction 

Electromagnetic calorimeters are often composed of inorganic scintillating crystals viewed by pho- 
todetectors. The energy resolution attainable depends primarily on the number of optical photons 
emitted by a scintillating crystal for a given energy deposit. Usually only a fraction of these pho- 
tons are collected by a photodetector. In the following we consider a photomultiplier (PM), where 
the collected optical photons are converted in photoelectrons with an efficiency characteristic of 
the photocathode and amplified up to the anode by the dynode system. 

Besides a possible non-linear crystal response [|TJ] , the anode charge pulse is therefore pro- 
portional to the energy deposited in the crystal by the primary particle and the fluctuations of this 
charge determine the energy resolution of the detector. The main contribution to this resolution 
comes from the statistics of the photoelectrons. In addition a small but not negligible contribution 
comes from the fluctuations in the PM gain and in particular from the gain of the first dynode. This 
resolution is worsened if the readout electronics integrates only part of the total charge delivered 
by the anode of the PM. This eventuality can occur when scintillating crystals with a long decay 
time are used in high-rate experiments where short integration times are needed. 

In the present paper we will consider the possibility of using BGO crystals and PM's for 
high rate experiments. In that case to limit the pile-up and the dead time, the PM pulses must be 
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integrated over a time interval shorter than the scintillation time of the BGO (Jscint = 300 ns). The 
deposited energy is then deduced from the maximum amplitude of the integrated pulse, measured 
by a peak-sensitive circuit. With this procedure only a fraction of the total charge is measured and a 
faster response is obtained at the cost of a worse resolution. This effect was experimentally studied 
for different integration times and the results were compared with a Monte Carlo (MC) simulation. 



2. Experimental setup 

The experimental setup is sketched in Fig. [j]. Photons with an energy of 662 keV from a 137 Cs 
radioactive source are detected bya2x2xl8 cm 3 BGO crystal, read at both ends by two EMI- 
9814B PMs. 
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Figure 1. Experimental setup. 



One of the two photomutipliers (PM-1) was used to trigger the acquisition of the pulses from 
the other photomultiplier (PM-2). In order to get rid of the noise, the trigger-signal from PM-1 was 
amplified and shaped with a gated biased amplifier (ORTEC-444) having an integration time of 
250 ns. The signal from PM-2 was processed by a filter amplifier (ORTEC-474) that has a variable 
gain and an integration 1 time (T, nf ) that can be set to 20, 50, 100, 200 and 500 ns. The output 
signal of that module was then acquired by a Lecroy WavePro 7300A digital oscilloscope, having 
a bandwidth of 300 MHz and a sampling rate of 250 MS/s. 

3. Equivalent circuit 

The ORTEC 474 integrating amplifier and its connection to the PM-2 anode can be represented by 
the equivalent circuit reported in Fig. @. When a given energy is released in the crystal at the time 
t = 0, the anode current 2 delivered by PM-2 is: 

/(f)=/oe~' /w x«(0 (3-1) 

'The differentiation control of the 474 module was set in the out position. This corresponds to a differentiation time 
of 0.2 ms which has a negligible effect on the output signals. 

2 The formulas reported in this section hold for a very large number of photoelectrons per pulse, so that the statistical 
fluctuations are negligible. 
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Figure 2. Equivalent circuit of the ORTEC-474 integrating amplifier. The two buffer amplifiers have an 
overall gain G while the integration time T,„ f is equal to RC. 



where T S cint is the decay time of the scintillator and u{t) is the unit step function. The total charge 
per pulse released by the anode of PM-2 is Q = Iq X sc i n1 . Using the standard Laplace transform 
method, the output signal V(t) of the circuit on Fig. ^ turns out to be: 

V (t) = _ e -th„„ ) (3 2 ) 

where /?,„ is the input resistance of the first buffer amplifier, G is the overall gam, Ti n t = RC is the 
integration time and a = %\ nt / T sc i„ t . For N pe — > °°, the maximum amplitude of the output signal is: 



A Npe ^=I R in a a l {{ - a) (3.3) 

which occurs at the time: 

In a 

T Npe ^o„ = %m -—^ (3.4) 
From Eq. the total charge (Q out ) of an output pulse is proportional to Q: 

Qout = I ^-dt = GI R in W = G R in Q (3.5) 

JO Rout 

4. Data taking 

About 15000 pulses were recorded for each of the five possible values %i nt = 20, 50, 100, 200 
and 500 ns. Each pulse was sampled every 4 ns during 10 jits and the resulting 2500 istantaneous 
amplitudes were acquired. A thousand pulses are shown in Fig. ^|a, while in Fig. 0b the average 
waveform is reported for the five T,„ f values. 

To check the reproducibility of the results, four data-sets (A, B, C and D) were acquired under 
different conditions for all values of T,„,. 

• A and B were taken under the same conditions (HV of PM-2 equal to 1 800 V and gain of the 
ORTEC-474 set to 10) but in different days in order to test the reproducibility of the results; 
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Figure 3. (a): A thousand of superimposed waveforms acquired with T,„ ( = 500 ns. (b): Waveforms 
obtained by averaging 15000 pulses: without integration (grey area not normalized) and with an integration 
time ii nt = 20, 50, 100, 200 and 500 ns respectively from the larger to the smaller signal. 



• C : HV of PM-2 equal to 1800 V and the ORTEC-474 gain set to 2; 

• D : HV of PM-2 equal to 1750 V and the ORTEC-474 gain set to 10. 

During the acquisition of the four data-sets the configuration on the trigger side of the set-up (PM-1 
and ORTEC 444) was kept fixed. The trigger level of the oscilloscope was set sufficiently low to 
accept all the pulses from the 662 keV photons. The electronics noise was measured by acquiring 
data with a random triggers. The width of these noise spectra being 10 times smaller than that of 
the source signal, the contribution of the electronics noise to the energy resolution was neglected. 

5. Data analysis 

The acquired waveforms were analysed off-line. For each value of Zj nt and for each pulse, the 
maximum amplitude (A), the peaking time (T) and the total charge (Q ou t) were evaluated. 

5.1 Total charge fluctuations 

When the total charge of each pulse is measured, the energy resolution of the BGO crystal is mainly 
determined by the statistical fluctuations of the total number (N pe ) of the photoelectrons and by the 
fluctuations of the gain of the 12 PM dynodes (gi,~-,gn)- Assuming a Poisson distribution for N pe 
and considering that in general gi > g2 = g3 = • •• = gn = g the energy resolution 3 is given by 

3 Throughout this paper x and o x indicate respectively the mean value and the r.m.s. of a Gaussian fit to the x 
distribution. 
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Figure 4. (a): A typical spectrum of the total pulse charge. The peak from 662 keV gammas is superim- 
posed to a background (grey area) which was measured without the 137 Cs source, (b): Relative total charge 
resolution as a function of T,„ r for the four sets (A, B, C and D) of data taking. The dashed line represents 
the expected behaviour for N pe = 130. The statistical errors are within the points. 
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where E is the measured energy and k = 12 is the number of dynodes. Since g is rather larger than 
1, Eq.|5J] was approximated as: 



(5.2) 



that is equivalent to take into account only the contribution of the fluctuations of the first dynode. 
In Fig. ^a a typical experimental spectrum of Q out is shown. Taking into account the HV of PM-2, 
the characteristics of the 9814B tube [J7|] and of the BeCu dynodes [Q] the average gain of the first 
dynode was assumed to be g l ~ 6. 
The energy resolution becomes: 



E Q out 




(5.3) 



In Fig. |]b the experimantal values of OQ out /Q ou t we reported for the five values of T;„ f and 
for the four sets of data taking. They show, as expected, a flat behaviour with respect to Tj nt . A 
constant fit to these data allows to determine from Eq. |53] the average number of photoelectrons 



N pe = 130^. 
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Figure 5. (a): Example of a spectrum of the maximum pulse amplitude A. The curve is a gaussian fit to the 
data, (b): Experimental maximum amplitude resolution (open points) as a function of T,„ ( . The crosses are 
the values predicted by Eq.'s 5.4 and 5.6 with N pe = 130 and with the measured values of T (Fig. ||). The 
statistical errors are within the points. The curve is the prediction of the MC (see Sec.Q). The dashed line is 
the asymptotic value of the curve, which is equal to Oq oui /Q ou ,. 



5.2 Maximum amplitude fluctuations 

When the measurement of the total charge Q takes too long, the energy deposited in the crystal can 
be inferred from the maximum amplitude A of the integrated signal. In Fig. |5]a a spectrum of A 
obtained in the present test with T, n , = 100 ns is shown. The fluctuations {a a) on A are obtained by 
a gaussian fit to the data. In Fig. ||b the experimental resolution Oa/A is reported 4 as a function of 
Ti nt . At large values of Xi„ t and in particular for x- mt » T sa ; lt , Oa/A tends to the value of OQ m jQ 0Ut , 
while at smaller values of T,„ ; the resolution worsen. To clarify the dependence of oA /A on T !nf , 
a naive Poissonian model based on an extension of Eq. |53] was adopted. According to this model 
the resolution is given by: 




(5.4) 



where n pe is the number of photoelectrons which contribute, for each event, to its maximum am- 
plitude A, i.e. those emitted before the peaking time T of that event. 

The average value n pe was approximated to the fraction F of the average total number of 
photoelectron emitted before the experimentally measured values of T: 

n pe = FN pe (5.5) 
F = (l-e~ r/w ) (5.6) 

In Fig. ^the measured dependence of T on T,, 1( is reported. 

In Fig. ^|b the predictions of Eq.'s and H3| are compared with the experimental results. 
While for Tj nt = 500 ns the agreement is good, at lower values of T,„ f the experimental resolution is 

4 The data from the four data sets A, B, C and D have been averaged. 
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Figure 6. Experimental average peaking time (T) for the five considerd values of T, n( . The error bars are 
within the points. The continuous curve is the prediction of the MC (see Sec.^) for N pe = 130. When 
N pe — > °° the continuous curve tends to the dashed curve represented by Eq. \5.4{ 

better than predicted with the naive model. To understand this discrepancy, a detailed Monte Carlo 
simulation of the experimental situation was performed. 

6. The Monte Carlo simulation 

The formulae reported in Sec. represent the response of an RC integrator excited by an exponen- 
tially decreasing current composed of a very large number of electrons, so that the charge quan- 
tization is washed out. In the situation we are considering, the average number of photoelectrons 
per pulse is relatively small so that the response of the device must be simulated with a MC. The 
input current I(t) is described as a sum of delta functions, each one corresponding to an incoming 
photoelectron. Then the amplitude of the integrated pulse at a time t turns out to be the sum of the 
contributions from all the photoelectrons emitted before that time 5 : 

n(t) 

V(t)=R in (^L)Y,G i e-( t - t ^"> (6.1) 

tint i= \ 

where q is the electron charge, n{t) is the number of photoelectrons emitted before the time t, ti 
(i = \ ,n) is the emission time of the z'-th electron (0 < < t) and G, is the PM gain for the z'-th 
electron. The probability distribution function of the is a decreasing exponential with a decay 
time equal to T sc int- For a fixed t,n{t) follows a Poisson distribution with a mean 

n(t) =N pe (l-e- ,lts ""<) (6.2) 

It is worthwile to note that Eq. ^TT] represents a single pulse that can be used to measure the 
energy, only if there is an effective pile-up of the contributions of many photoelectrons belonging 

5 It can be shown that Eq. ^H] — !> Eq. when N pe — > °°. 
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Figure 7. (a) A typical pulse generated by MC with N pe = 130 and T,„ r = 500 ns. Each sharp discontinuity 
in the pulse shape corresponds to the arrival of a photoelectron. (b) Example of experimental pulse recorded 
with the same T, nr . 

to the same detected particle. This occurs if the integration time %\ M is much larger than the average 
time interval between two consecutive photoelectrons z sc i nt / 'N pe )\ 



while if K < 1 the energy released in the scintillator gives rise only to a series of single photoelec- 
tron pulses. In the present experiment K ranges from 8.7 (at T,„ f = 20 ns) to 217 (at Xj nt = 500 ns). 

In the MC simulation all the aforementioned effects were taken into account. The MC was 
run with 10 < N pe < 10 4 and for %- mt ranging from 10 ns to 1 /is. For each of the ~ 100 pairs of 
values {N pe , T,„;) about 10000 pulses were generated. For each simulated pulse the MC calculates 
its amplitude V(t) every ns during an interval of 2 /is. The maximum amplitude A and the time 
T at which this maximum occurs were recorded for each pulse and the relative fluctuations Ga/A 
were determined. 

6.1 Comparison with experimental data 

For a comparison with the experimental data, the MC was run with N pe = 130. The Poisson 
fluctuations of the gain g of the first dynode, with a mean value g = 6, were also taken into account. 

In Fig. |7J a typical pulse generated by MC with Tj nt = 500 ns is compared with an experimental 
pulse recorded with the same integration time. The agreement between the two shapes is quite 
good. 

In Fig. |5]b the dependence of Oa /A and in Fig. || the mean value of the peaking time distribution 
(T) on Tim, calculated with the MC, are compared with the experimental points. In both cases the 
agreement is quite good. This confirms that in the present experimental conditions the experimental 
resolution is better than predicted by the naive model based on a Poissonian statistcs. These checks 
give confidence in the MC simulation and allow to use it to predict, in the most general experimental 



K = N i 




> 1 



(6.3) 



situation, which is the energy resolution attainable with an integrator followed by a peak-sensitive 
electronics. 

6.2 Energy resolution in the general case 

The resolutions ox/A, calculated as a function of a = tint/tstint and for different values of N pe , 
are reported in Fig. |8] which is therefore a general utility to evaluate the resolution attainable with 
a scintillator having a decay time T iC! „,, read by a photodetector followed by an integrator and a 
peak-sensitive electronics. 

To present these results in a general form the fluctuations on the photodetector gain have 
not been included because they depend on the particular type used. For a PM the effect of these 
fluctuations can be taken into account by multiplying the values of ox/A read on Fig. || by the 
corrective factor of Eq. pTTj . From Fig. ^| it appears that to perform an integration with T,, ? , < T sc i nt , 
at least 10 photoelectrons are needed. 

As already pointed out in the naive model, only the n pe photoelectrons emitted before the 
peaking time T contribute, for each event, to the the maximum amplitude A. Assuming a Poisson 
distribution for n{T) the relative fluctuations on A is given by: 
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Figure 8. Look-up figure, calculated with the MC simulation, which allows to evaluate the maximum 
amplitude resolution for different values of the integration time, and for different values of N pe . Only the 
region where K > 5 (see Eq.|b.3D is shown. 
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Figure 9. (a) Relative maximium amplitude resolution as a function of the inverse square-root of the 
average number of photoelectrons which contribute to the maximum amplitude of the signal. For the first 
and last points of the upper three series the value of a = T,„ r / T SC mt is reported, (b) Average peaking time 
divided by the peaking time at N pe — > °° (Eq. 3.4) as a function of the actual number of photoelectrons in 
each pulse, for different values of a. Only the region where K > 5 (see Eq.|6.3D is shown. 



Contrarily to the experimental situation, in the MC simulation n pe is a known quantity for each 
event, so that the simple model represented by Eq. SA can be tested. In Fig. ^a cr^/A, calculated 
with the MC, is reported as a function of 1/ ^/n pe , for different values of N pe and a. It appears that 
for Np e > 500 Eq. 6.4 is satisfied for any of the considered a values so that the statistics of n pe is 



Poissonian and the naive model is valid. For smaller values of N pe and in the range where Eq. 6.3 



is satisfied, the resolution is better than predicted by Eq. |T4] so that in that region the statistics is 
sub-Poissonian. The reason of this behaviour is clear from the curves reported in Fig. 0b: for a fixed 
integration time a positive (negative) variation of N pe with respect of its average value N pe results 
in a negative (positive) variation of the corresponding peaking time, which partially compensates 
the variation on N pe . This anticorrelation between N pe and T is responsible for the sub-Poissonian 
fluctuations at the lower values of N pe . 

The results reported in Fig. || and Fig. |9| are valid for any integration time and for any decay 
time of the scintillating light, when the readout electronics measures the maximum amplitude of 
the integrated pulse. 



7. Conclusions 

The possibility of using a scintillating crystal with a slow decay time (like BGO) for an electro- 
magnetic calorimeter in a high-rate experiment was investigated. In these experimental conditions 
a fast measurement of the energy deposited in the crystal is needed. This can be obtained, at the 
cost of a lower energy resolution, by integrating the output signal of the photodetector over a short 
time and by acquiring the maximum amplitude of the integrated signal. 
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An experimental test and a Monte Carlo simulation show that the energy resolution comes 
from the statistics of the number of photoelectrons emitted before the peaking time of the 
integrated pulse. While for a large number of photoelectrons the statistics follows a Poisson 
distribution, at a lower number of photoelectrons the statistics becomes sub-Poissonian due to an 
anticorrelation between the fluctuations of the number of photoelectrons per pulse and the peaking 
time of that pulse. The results are reported in a general form which allows to evaluate the 
contribution of the photoelectron statistics to the resolution of a calorimeter equipped with a 
scintillating crystal read by a photomultiplier, followed by an integrator and a peak-sensitive 
electronics. 
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